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DISCUSSION 


deduce the gauge-discharge relation river seems over 
optimistic. The following phenomena, that seem well established, 
illustrate the reasons for the writer’s belief: 

laboratory flume with bed natural regime-type 
river sand, constant discharge, and various bed-sand injection rates, observa- 
tions after steady regime depths had been attained showed that: 


Rugosity was greatly altered changes injection rate within the 
practical large-river range; 

Bed-wave full amplitude was the order one-fifth water depth; 

After moderate displacement from regime, 30-ft flume, carrying 
about one-fifth per sec, requires some settle down regime; 

Bed-wave conformation changes markedly during attainment regime, 
and bed waves have greatly enhanced amplitude when degradation occurring. 


large regime-type daily sonic observations bed waves 
during month including high flood peak showed that: 


Bed-wave length closely related discharge; 

Bed-wave full amplitude 40% water depth high discharge. 

Mean size bed material reach between diversions very closely 
correlated with estimated formative depth. (Samples were taken daily 
various verticals which formative depth did not vary appreciably with time.) 

Discharge division between two branches found agree reasonably 
with the regime theory slope formula terms bed-sediment size. 


the same river was found that distribution sand sizes within 
bed samples follow simple universal law. 

river models has been found that successful models have 
width and depth scales related consonance with regime theory derived from 
the study self-formed, sediment-bearing canals controlled discharges. 


terms the preceding facts the following phenomena might expected 
rivers: (a) Rugosity dependent rate change discharge, that is, 
hydrograph form; rugosity dependent bed-load discharge intensity, 
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which can affected conditions upstream; (c) rugosity affected meander 
pattern; (d) some average discharge, conditions explicable terms 
regime theory for steady discharge; (e) Eq. which pertains noncorrugated 
boundaries, will inapplicable except after severe distortion its 
and (f) for the conditions which regime theory applies, only one parameter 
should needed define the bed-sediment size. 


Therefore, the points Fig. should expected scatter about mean 
line representing regime conditions. The coordinates are nondimensional 
numbers containing the principal pertinent variables with the exception 
bed-load discharge intensity. The line cannot expected have simple 
equation since has correct for the misuse Eq. the two values 
are chosen arbitrary manner; and the use any kind D-term presum- 
ably hides the settlement velocity that does not follow simple law the 
size range river sands. 

Fig. considerable value the writer. shows that, choosing 
suitable river sites and analyzing data scientifically, regime-type relation 
demonstrated, just such relations were shown the much simpler data 
controlled canals. Considering the scatter points canal analyses when 
regime theory was being established, the scatter But, 
instead using the points Fig. establish the existence regime 
relation, the smoothing curve used indicate where the point for given 
stage given river likely be, then the scatter very bad; for example, 
where the curve indicates value the correct answer appears 
anywhere between about and depending host circumstances 
which there measure. 

The writer emphasizes that his criticism application; holds 
Dr. Einstein’s basic work, from which the present paper springs, the highest 
regard and has used practical applications. 


applying the Manning formula natural stream, engineers well 
engineering students often ignore two important facts: 


natural streams the flow rarely uniform because the effect the 
inevitable change discharge due the varying shape and size the channel 
cross section and the unavoidable irregularities channel alinement; 

natural streams the roughness value not constant, varying with 
the discharge and the depth water the channel. This can best illus- 
trated typical case, such shown Fig. indicating how the roughness 
varies with the depth water the channel and with the corresponding varia- 
tion discharge. The n-value high high stages because rough and 
grassy banks and increases rapidly near the bottom where debris has accumu- 
lated. 


Engineers and students usually learn the Manning formula from college 
hydraulies, which the theoretical cases are treated (for which cases the formula 


Prof., Hydr. Eng., Univ. Urbana, 
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designed only for uniform flow and single-valued generally assigned), 
often happens that practical problems they apply the Manning formula 
indiscriminately natural streams the same way theoretical artificial 
channels. The condition uniform flow and constant roughness uncon- 
sciously assumed without apparent reason. The writers are glad that the paper 
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points out that the hydraulic computations for roughness natural channels, 
the assumption uniform flow may made with negligible error. practical 
problems the variation with evident that there need develop 
new and more rational procedure for estimating the varying roughness 
natural channels. 

seems that the procedure hydraulic computation described the paper 
may simplified with practical advantages that the most familiar form 
the Manning formula still retained use and that several terms unfamiliar 
practicing engineers are eliminated. 


, 
Relationship Between and for Given Values the paper 


the definition and and the equation for are given follows: 


and 


58.68 


Eliminating and from the equations and simplifying, 
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From Eq. 17, 
D yp’ 


Now, assuming various values R’/R, Eqs. and will give the rela- 


the relationship between V/V’’; and shown the curve Fig. the rela- 


tionships between and and between and y’, respectively. 


tionship between and can obtained indirectly. For example, 


assume 0.60, then, from Eqs. and 24, and 


Take 5.60. Then, from Fig. 11.5. The corre- 


sponding values 1.79 and 2.00. The writers have com- 


puted number pairs the values the item and for various 


values R’/R. The results are plotted the form curves shown Fig. 
Equation for Roughness Factor the Manning formula 


(3) 
and Eq. gives 
which 
(20) 


The roughness factor may then expressed eliminating V’;, and 
from 20, and 25: 


Procedure Hydraulic means the curves Fig. and 
Eq. the practical procedure for computing the rating curve may performed 
follows: 
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Values 


CURVES 


0.2 
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Using the computed values and find the value from 
RS’ 
Fig. 
Compute Eq. 26, using the given and the value obtained 
step and 
Compute the discharge the Manning formula. 
6.7 3.81 13.24 0.688 0.260 0.0286 1,010 2,860 
8.1 5.00 14.50 0.524 0.366 0.0229 5,850 
9.0 5.85 15.25 0.448 0.513 0.0178 


Using the authors’ example for illustration, the procedure is: (1) Assume 

3.81, (2) 13.24, and 0.688; (2) from Fig. 0.26; 
(3) Eq. 26, 0.0286; and (4) using the Manning formula, 2,860 


per sec. The complete computation shown Table The corre- 
sponding points involved step are indicated dots Fig. 
clear that the idea presented this paper developed into procedure 
that practicable and workable. However, the accuracy the procedure 
depends many factors. The two most important ones are: 


(a) Satisfaction with the use describe the grain roughness. 
cording the theory statistics the distribution observed variate may 
described number parameters. engineering works, parameters 
greater than the third order (skewness) are usually insignificant, and hence 
rarely employed. However, the other hand, single parameter too in- 


TABLE VARIATION COMPUTATION 


DIscHARGE CuBIC 


centage 
vari- 

Maxi Mini Maxi Mini Maxi Mini 


mum mum mum mum mum mum 


1.00 8.25 5.40 6.70 4,020 1,415 2,860 0.044 0.023 0.029 
2.00 9.60 7.00 8.10 8,160 5,030 6,200 0.025 0.019 0.023 26.1 
3.00 10.00 8.10 9.00 11,600 8,120 9,630 0.020 0.017 0.018 


definite describe distribution. seems that the distribution grain size 
irregular that cannot described completely any single value 
representative diameter unless the pattern distribution the same type 
and shape for all possible cases. 


— 


Discharge, Thousands Cubic Feet per Second 
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(b) The accuracy the theory upon which the curve Fig. based. 
The authors have obtained the working relationship between the values 
from actual measurements for number rivers. This relation- 
ship represented smooth curve fitted the observed data. However, 
appears that the plotted observed points scatter within band certain 


Depth, Feet 


width rather than along narrow determinable path that could well repre- 
sented aline. interesting see what the effect this scatter the 
computed roughness factor Table gives computation which the ex- 
treme values taken from the band the plotted points Fig. are used. 
The range variation amounts high 72.4%. The resulting rating 
curves are shown Fig. 


ASCE.—With rare exception, the historical 
problem hydraulic roughness has been concerned with the development 
empirical expressions for determining fluid resistance. However, the case 
which one type roughness may considered superimposed over 
another type, empirical formulas expressing the friction loss the total 
roughness either not exist, the authors’ problem, the formulas in- 
troduce such vague characteristics “equivalent sand Con- 
sequently, the authors remark, the determination friction losses this 
case necessarily based entirely upon judgment. Present-day knowledge 
the fluid flow should certainly warrant more precise statements 
tortuous river with considerable vegetation possesses Manning 


coefficient The rational approach proposed the authors the 


Fulbright Fellow, Imperial Science and Technology, London, England. 
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problem river channel roughness thus permits the engineer calculate the 
Manning coefficient for this kind roughness. 

The writer, conjunction with Mr. Einstein, conducted experimental 
investigation concerned with the resistance exerted boundary composed 
two more types friction-producing elements acting over the same bed 
Results obtained from this investigation have led the conclusion 
that other types hydraulic problems also may solved utilizing the 
concept given the authors. For example, possible calculate the 
Chezy coefficient for particular flow for the case which small sills are at- 
tached, widely separated, flow boundary. this case, one assumes that 
the total resistance the flow equal the sum the resistances exerted 
the sills, and the boundary area between the sills. Consequently, 
expressions are available which give the resistance exerted each type 
roughness, the total resistance can obtained directly, and from the latter 
the Chezy coefficient computed. 

One might consider the authors’ parameter, ratio weight shear. 
Multiplying the numerator and the denominator Eq. seen that 
the former becomes the submerged weight layer height, and that the 
latter becomes the grain shear. Suppose that flow increased; then the 
rate bed load transport increases and, accordingly, decreases. For 
decrease, the denominator Eq. must increase, which compatible 
with the assumption regarding increased flow. 

hoped that river measurements taken others the future can 
calculated terms the ordinates Fig. and included with the authors’ 
data. Meanwhile, employing the curve which the authors regard tenta- 
tive and utilizing their concept resolution shears, the problem pre- 
dicting river channel friction losses can solved. 


flow” formulas not permit 
accurate computation discharges for varying stages. There are two ways 
adjusting such formulas produce accurate results. One adjust the 
roughness factor with changing stage and the other change the river cross 
section. Since rules have existed for making such adjustments, the pro- 
posal this paper will very helpful. 

The writer concurs with the authors (discussion preceding Eqs. 18) that, for 
practical purposes, the Manning-Strickler equation (Eq. 6).is more convenient 
use than Eq. Eq. because solutions are direct. However, 
desired use Eq. practical application, the data used develop the 
relationship (as the curve Fig. should computed use Eq. and not 
Eqs. and was done the authors. 

The writer suggests the following procedure evaluate the necessary re- 
lationships for the Manning equation: Assume that the Manning roughness 

Am. Geophysical Union, August, 1950, pp. 603-610. 

Friction Open Channel Flow,” Banks, dissertation presented the 
Univ. California, Berkeley, Calif., June, 1951, partial fulfilment the requirements for the 


degree Doctor Philosophy. 
Hydr. Engr., Garrison (N. Dak.) District, Corps Engrs., Dept. the Army. 
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factor composed two parts: 


which 


roughness factor depending channel bed grain diameter and the 
roughness factor depending channel bed irregularities. 

The Manning equation expressed Eq. and the Manning-Strickler 
equation 


which identical with Eq.6. The expression for given Eq. properly 
expression for n’, the grain roughness, and the roughness factor Eq. may 
and that 


3/2 
(30a) 
or 


reasoning similar that set forth the authors the discussion 
preceding Eq. 17, assumed function the sediment transport 
determine the relationship between and the data from actual river 
flow measurements were used—namely, all data which Fig. based sup- 


0.070 


Values 


Fra. 


plemented those for the Missouri River, near Williston, Dak. The 
following line computations was used: from 28, from Eq. with 
and from flow measurements, from Eq. 27, from Eq. 30, and 
from Eq. adjustments were made for friction losses banks. 


0.030 
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The results are indicated Fig. with the ordinate and the 
abscissa, indicating the increase sediment transport from left right the 
same manner that proposed Carl Brown,” ASCE. Fig. the 
roughness factor that reflects the channel irregularities n’’ decreases with in- 
creasing flow and has minimum value approaching zero. 

The procedure for practical application could the following: Compute 
from Eq. 28; assume R’; determine from Eq. 17; read off from Fig. 
and calculate from Eq. 27, from Eq. 30, and from Eq. repeated 
use, advisable draw curve with ordinate and abscissa. 

Acknowledgment.—The co-operation Mr. Einstein providing the flow 
data used constructing Fig. and Pratt, ASCE (Chief, Reports 
and Hydraulics Branch Garrison District the Corps Engineers), for 
suggestions and for permitting the use the data for the Missouri River, near 
Williston, gratefully acknowledged. 


ASCE.—The physical approach the subject 
channel roughness and its relationship sand transport sound long 
the authors restrict two-dimensional flow; but over-simplification 
assume that all energy losses except those due grain roughness were caused 
shear force resulting from irregularities the periphery. laboratory 
investigation energy losses one-way (river) and two-way (estuary) flow, 
the Hydraulics Research Station the Department Scientific and In- 
dustrial Research Wallingford, Berkshire, England, the energy losses have 
been divided into three categories: (1) Those due textural roughness; (2) 
those due ripple roughness; and (3) those caused form drag resulting from 
the major irregularities the banks, bed, islands, and sandbanks—all which 
induce curvature flow and tend produce meanders. The authors cannot 
justifiably group these losses with irregularities the bed due sand ripples, 
bars, and sand dunes because, whereas the bed pattern and bed movements 
adjust themselves fairly rapidly the discharge varies, this does not occur 
the case major irregularities such promontories, outcrops, gradually 
changing banks, islands, and sandbanks, which have entirely different time 
from that bed pattern adjustment the movement sand grains. 

regards the size grains exposed the bed—the writer cannot 
agree that the size remains constant and that the 35:65 percentage distribution 
material remains the same all discharges. India, the writer found ma- 
terial suspension, material exposed the bed, and also bed levels vary 
considerably with different discharges. collect samples bed material 
various discharges, the other hand, would make the work very laborious. 

also appears that the authors have assumed not merely the same drop 
between the upstream and downstream observation points but also the same 
uniform slope for all discharges. Actually, however (as the authors themselves 
state), the curvature flow and, hence, the distance the water travels, de- 
creases the discharge increases. other words, the slope increases; and 
“Engineering edited Rouse, John Willey Sons, Inc., New York, Y., 1950, 


Director, Hydraulics Research Station, Dept. Scientific Industrial Research, Wallingford, Berk- 
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many rivers the slope far from uniform, lower stages consisting series 
pools and rapids. 

Rather than rational approach, the writer considers that linking the flow 
with grain size which assumed unaffected discharge, advanced 
this paper, extreme empiricism, which variables are assumed con- 
stants. 

The next question concerns the degree accuracy attained. Crump 
(senior engineer the writer’s staff) has analyzed the data. Be- 
ginning with the first item Table (Missouri River, Pierre, Dak.), 
found trial and error the results shown Table The labor involved 


R’ R” R’ R” R” /R’ Q Qc 


Stage 
puted sured puted sured puted sured puted sured (%) (%) 


obtaining these results was great that Mr. Crump adopted simpler method 
for obtaining the probable error, which had the added advantage that repre- 
sented average performance the several channels the same manner does 
the standard line Fig. which the authors for solving specific 
problems. Disregarding results for the Salinas River Paso Robles, the 
scatter results has standard deviation 20% (very approximately) 
(1.68 
are known, and the values and for various values the gage have been 
extablished survey, the authors’ method calibrating the gage for discharge 
assume successive values the component R’. For each value 
assumed, the velocity determined from Eq. and for the known value 


Fig. given site where the constants K,, and 


2 2 


reading the value from the curve Fig. the value 

2 2 

Following this method Mr. Crump found that the average value 
was 5.80 for low discharges, 2.00 for intermediate discharges, and 0.66 for high 
and, accepting the foregoing average values R’’/R typical, 


may range from 


= 
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The effect these standard deviations values shown the 
ratios and difficulty expressing the effect terms rela- 


tive discharges will apparent from Fig. which indicates the relationship 


TABLE THE EFFECT STANDARD ERROR 


(1) (2) (3) (4) (5) (6) (7) (8) (9) 
5.80 6.80 1.562 0.694 10.062 5.025 1.480 0.739 
Medium....... 2.00 3.00 1.562 0.694 4.124 2.388 1.375 0.796 
0.658 1.658 1.562 0.694 2.028 1.456 1.223 0.877 
Explanation: 


the headings, subscript denotes values given the standard curve, Fig. and subscripts 0.8 and 
1.2 denote values that are, respectively, 20% less than, and 20% greater than, those given Fig. 


between and for three values 0.8, 1.0, and 1.2, re- 
spectively. The curves are drawn for the case medium flow, for which the 
0.796. Being based common value the component R’, all such points 
have the same value and therefore lie horizontal required 


Values 


Values 


know the ratio the ordinates points and which lie the same 
vertical and therefore have the same value (gage reading) and (sec- 
tional area) that the ratio the discharges the same the ratio 


| 
j 
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the velocities know these ratios necessary know the 
value the index the assumed relationship Fortunately, 
Table the authors have given two values and seven the nine sites 


taken into consideration. With the value given 
2 2 
1.27 7.72 0.84 


Average 


Computed values and for the seven cases considered are given 
Table For the authors’ results, Table these computations establish the 
average relationship, 

(31) 


Comparison with Manning’s formula— 


(32) 
—requires that 
0.74 


Eq. yields 0.74 against 0.755 found analysis. 

Reverting Fig. 10, with now possible find the ratio the 
velocity ordinates (which also the ratio relative discharges) the three 
points and having the same sectional area. moving along the lower 
curve from points (Fig. 10) the abscissa reduced the ratio 1.375:1, 
that the ordinate reduced the ratio 1—that is, the ratio 
1:0.639. This also the ratio the discharges, the percentage 
error discharge being 36.1%. Similarly moving from point 
point along the upper curve, the ordinate increased the ratio 
is, the ratio 1:1.381; that 1.381 giving error 
38.1%. This gives the result for medium stage flow. Results for the 
three stages flow are shown Table 


| 
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the foregoing analysis, Mr. Crump has shown that the errors resulting 
from following the authors’ method are far excess the 10% error which 
experienced river engineers India were able work assessing the dis- 
charges selected short, straight reaches, and the claim (see sentence preced- 


FROM 


= 0 80 = 1.20 
Stage flow 
(1) (2) (3) (4) (5) (6) (7) 
1.480 0.575 0.739 1.532 +53.2 
Medium..... 1.375 0.639 0.796 1.381 +38.1 
1.223 0.753 —24.7 0.877 1.203 +20.3 


ing that obvious advantage this method that 
the results not depend upon individual judgment for selection the appro- 
priate shown unjustified; because the errors due the 
author’s assumptions far exceed the errors anticipated assessing the 
appropriate value for selected reach. 


yar 
é 


